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The phase structure function accumulated by two parallel rays after transmission through a layer 
of turbulent air is best known by a proportionality to the 5/3rd power of the lateral distance in the 
aperture, derived from an isotropic Kolmogorov spectrum of the refractive index. For a von-Karman 
spectrum of the refractive index, a dependence involving a modified Bessel function of the ratio of 
the distance over the outer scale is also known. 

A further standard proposition is a proportionality to the path length through the atmospheric 
layer. The manuscript modifies this factor through a refined calculation of an integral representation 
of the structure function. The correction establishes a sub-linearity as the lateral distance grows in 
proportion to the layer thickness; it is more important for large than for small outer scales. 
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I. SCOPE 

_i ^ A. Planar Phase Screen Terminology 
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A phase screen characterizes the distribution of optical path lengths in the aperture of a receiver depending on the 
inhomogeneity of the refractive index along the path through the atmosphere up to the emitter [ll[. The phase tpir) 
in the receiver's pupil plane is a function of the two-dimensional spatial vector r. The phase is the line integral of the 
^ • product of the vacuum wave number fc = 27r/A by the index of refraction n along the path through the atmosphere 
at the wavelength A of the electromagnetic spectrum. 

The seeing quality, figures of merit like the Strehl ratio, and the requirements on adaptive optics systems aiming to 
straighten these wave fronts refer to these phase screens. The current work is a contribution to an accurate description 
: of the phases, the measured quantities, in terms of parameters that describe the turbulence (structure constant and 
outer scale length) and in terms of geometric parameters which determine the sampling (optical path length of the 
lO ! rays and distance between these). 

'/"J ' In the remainder of this section, the manuscript recalls how a model of isotropic turbulence is set up, then reviews 
\ in Section [IT] how this leads through a procedure of integration of the refractive index along two parallel rays to a 
1 standard formula of the phase structure function of von-Karman turbulence. This is redone in Section [IIII on a more 
' accurate basis that keeps the layer thickness parameter consistent throughout the computation, with the main result 
that a finite layer thickness generates a smaller phase structure function than with the standard formula. 



B. Three-Dimensional Turbulence Spectrum 

This section summarizes the basic theory which links the structure functions of three-dimensional indices of refrac- 
tions to phase covariances measured in pupil planes. 

The structure function of refractive indices n is defined as the expectation value of squared differences if measured 
at two points separated by a 3-dimensional vector Ar, 

= ([n(r)-n(r + Ar)]2). (1) 

Assuming that the expectation value of the mean does not depend on r, binomial expansion rephrases this as a 
background value minus two times a correlation function, 

X'„(Ar) =2(71^) -2(n(r)n(r + Ar)). (2) 
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Its Fourier Transform as a function of spatial frequencies / is 

Vn{Ar) = 2 j $„(/)[! - cos(27rf • Ar)]dV- (3) 

The imaginary part oc isin(27rf • Ar) vanishes as we assume that the structure function has even parity: I?„(Ar) — 
!?„(— Ar). Other notations emerge if the factor 27r is absorbed into the spatial frequency or if the normalization of 
the Fourier Transform is chosen differently. 

We focus on the von-Karman model of an isotropic power spectrum with zero inner scale, a prefactor c„C^ — which 
is split into a constant c„ and a structure constant — , an outer scale wave number /o representing the inverse of 
the outer scale length, and a power spectral index 7, 

= c„C^(/2 + /„2)-(-+3)/2. (4) 

The Kolmogorov Model is obtained by setting 7 = 2/3 and setting the outer scale to infinity, i.e., /o = 0. This $ is 
inserted in Polar coordinates are introduced in wave number space with Jacobian determinant sm.(9 f), and 
f • Ar — frcosOf. Integration over the azimuth in these coordinates yields a factor 27r, and the integration over the 
colatitude 0f introduces a shape factor with a spherical Bessel function jo{x) = sm{x)/x, 

I?„(Ar) = iTTCnClJ ^ /2)(^+3)/2 >< [1 ^ cos(2^/r COS sin(0/)d0/ (5) 
f°° f ^ df 

= ^^c^Glj^ (^,/^2)(,+3)/2 [^-Jo(27r/r)]. (6) 
We require that 2?„ does not depend on c„ in the Kolmogorov limit: 



and this fixes the constant c„ as [la. 111 



in particular c„ « 0.009693 if 7 = 2/3. 



limP„-^C^^ (7) 

JO— s-u 



r[(3 + 7)/2] 



27r3/2+7r(-7/2)' 



II. LINE-OF-SIGHT PATH INTEGRALS 

A. Phase Structure Function 

We review the integration of isotropic structure functions I?„ along the line of sight [H, p. 293] [1, (Cl)][ll!, (7a)]. 
The distance to the reference point of the structure function is split into a horizontal vector component b (baseline) 
and a scalar component h. The phase Lp of the electromagnetic wave accumulated after transmission through the 
turbulent layer is the product of the optical path length by the vacuum wave number k, and the optical path length 
is the path integral over the product of geometric path length and refractive index n along the atmospheric height h: 

f-K j sin a 

V3(b) = fc / n(b,/i)d/i (9) 

Jo 

The scale height K is the vertical thickness of the turbulent layer. The air mass is 1/sina as a function of the star 
elevation angle a above the horizon. The structure function of the phase is defined as 

= (1^(0) -^(b)|2). (10) 
The expectation value of the square is expanded with the binomial theorem. 



i^/sin(a) f^Kj sin(a) 

kn{0,h)dh— / kn{b^h)dh 
^0 

pK/ sin(a) p^/ sin(a) p^/ sin(a) f^/ sin(a) 

2k^( I n{0,h)dh / n{0,h')dh'- / n{0,h)dh / n{b,h')dh' 

Jo Jo Jo 



(11) 

(12) 
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We switch to relative and mean coordinates, A/i = h — h' and h = {h + h')/2, with inverse mapping h = h + Ah/2, 
h' = h-Ah/2, 



2k' 



K/(2sin{a)) _ r2h 

dh / dAhn{0,h + Ah/2)n{0,h- Ah/2)) 

J-2h 
K/sin{a) _ r2lK / sin{a) -h] 

dh / dAhn{0,h + Ah/2)n{0,h- Ah/2)) 

K/(2sin{a)) J 2[h-K / am(a)] 

dh / dAhn{0,h + Ah/2)n{b,h- Ah/2)) 

J-2h 
K/sin{a) j.2lK / sin{a) ~h] 

dh / dAhn{0,h + Ah/2)n{b,h- Ah/2)) 

K/{2sin{a)) J 2[h- K / sin(a)] 



(13) 



In the first two terms of the right hand side, the distance between the locations in the arguments of the refractive 
indices is Aft., because these are located on the same ray. In the remaining two terms of the right hand side, the distance 
is between points located on two different rays, horizontally translated by an amount 6, rays slanted by the altitude 
angle a towards the horizontal. The location (0, h + Ah/2) on the first ray translates into Cartesian coordinates {h + 
Aft/2) (— cos A cos a, sin A cos a, sin a) in a suitable topocentric system with star azimuth A. The location (6, h — Ah/2) 
on the second ray has Cartesian coordinates (ft — Aft/2)(— cos Acos a, sin Acos a, sina) + 6(— cos Af,, sinAfc, 0) where 
Ah is the baseline azimuth [14i|. The square of the distance between these two locations is + (Aft)^ — 2bAhcos9, 
where cos 9 = cosacos(j4;, — A) measures the angle 6 between baseline and pointing. 




FIG. 1: Sketch of geometric parameters. The sections of the two parallel rays are the bold lines, a is the angle between the 
rays and the horizontal. The vertical distance between the top and bottom of the layer is K; so the geometric length of each 
ray is K/ sin a. 



In Figure (TJ is the angle between b and what is called the delay D in astronomical interferometry, 

D ~ bcosO, 

and the third leg in this rectilinear triangle is the projected baseline 

P = 6 sine. 

Reverse use of 

Vn = 2(n(0, ft)2) - 2(n(0, h)n{b, ft)) 
and insertion of the aforementioned distance b^ + (Aft)^ — 26Aftcos6' ^ + (D — Aft)^ yields 



2h 



' / dh dAftX>„(v/p2 + {D~ Aft)2) 

Jo J-2h 



(14) 
(15) 
(16) 



pK/{2sin{a)) _ p2h _ _ pK/(2sin(a)) 

2P(/ dh I dAhn^{0,h)) -P dh I dAhV^{Ah) 

Jo J~2h Jo J-2h 

pK/sm{a) p2[K/ s\n(a)-h] pK/sm{a) p2[K / sin(a) -h] 

+2P(/ dh I dAhn^{0,h)) ~P dh dAhVniAh) 

JK/(2sin(a)) J2[h-K/ sin(a)] Ji^/(2 sin(a)) J2[h-K/ sin(a)] 

pK/{2sin{a)) _ n2h _ _ n K / {2 sin(a)) _ f2h 

-2P(/ dh i dAhn'^{Q,h)) +P dh 

Jo J-2h Jo J -2h 

l-K/s\n(a) _ i-2[K / s\n(a)-h] _ _ pK/sinia) _ <.2[A:/ sin(a)-S.] 

-2P{ dh dAhn^{0,h)) +P dh dAhV^i^^P^ + [D - AHM 

JK/{2sin{a)) J2[h~K/ sin{a)] J_fC/(2 sin(a)) J2[/i-if/ sin(a)] 
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The four expectation values of the squared mean values cancel, 

r.A:/(2sin(a)) /.2h 



= -P / dh dAhVn{Ah) 

Jo J-2h 



pK/sin(a) _ p2[K / sin{a) - h] 
/ dh dAhVn{Ah) 

Jif/(2sin(a)) J 2[h- K / sin{a)] 

rK/{2sin(a)) /.2h 



-k 



/ dh dAhV,,{^ P"^ + {D - Ahf) 

Jo J -2h 



J -2h 

K/sin{a) f2lK/ sin{a)-h] 



+ / dh dAhVni^/P^ + {D~- Ahf). (18) 

JK/(2am(a)) J 2[/i-K/ sin a] 

In the second and fourth term substitute h for Q = K / sin(a) ~ hio observe that their contributions double the values 
of the first and third, 

f.if/(2sin(a)) ^2E 



= -2k^ I dh j dAhVn{Ah) 

Jq J -2h 

p K / {2 sin{a}) _ <.2/i 

■ / dh dAhV„{y/P^ + {D ~ Ahy) 

Jo J ~2h 



^K/{2ain{a)) /.2h 



= 2P / dh dAh VnWP^ + {D- A/i)2) - Vn{Ah) . (19) 

Jo J ~2h L 

Insertion of the Wiener spectrum © and substitutions Ah — Py and h = Pz/2 of the integration variables transform 
this into 



= 167rc„C2pP / ,2l+3)/2 / / jo{27rfPy)-jo{27TfP^l + {h cos 9 / P ~ yfllQ) 

Jo [J +Jo) Jo J-2h/P L J 

i-oo f'^df i-K/iPsina) i-z 

= SnCnClPP^ + f2^iy+3)/2 I dzj dy[joi2nfPy)-joi2nfP,/TTiD/P^)\- (21) 



B. Review of the Standard Theory 

The standard proceeding assumes that the contributions of the two spherical Bessel functions in (pil)) cancel well for 
large y. In consequence, the limits of the innermost integral are extended to h/P ^ (X), emitting a Bessel Function 
of order zero [i,[ll[2ll: 



/oo poo 
dyjo{27TfPy) - / dyjo{27rfP^l + [bcose/P-y]^) 
-OO J —OO 



-2JP-2JP'"^'''^''^- (22) 

This reduces the integration over h to a, simple factor d, ^ 

The remaining frequency integrals are known in the literature [lol 6.565.4] the frequency-to-real-space matching factor 
([8|) is also inserted, and the structure function for the von Karman spectrum is condensed into the well-known 0, 0] 



sin a " r(-7/2)(7rP/o)i+^ 



r (^^) -2(7rP/o)(i+^)/^if(i+^)/2(2^P/o) 



(24) 



where K (.) are modified Bessel Functions. [1, 9.6] The factors that depend on the projected baseline (phase) measured 
in units of the outer scale, 

Po = 27r/oP, (25) 
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are the topic further below. The Kolmogorov hmit of infinite outer scale is 



where 



5(7) ^ V^ ^^ " 2.9143808 



(26) 



(27) 



at 7 = 2/3. In overview, observations away from the zenith multiply the sale height with the air mass, and employ 
an effective distance between the two beams (which scales oc 5/3) equal to the projected baseline. 

III. CORRECTION FOR FINITE LAYER THICKNESS 

A. Aim 

Quantitative understanding of effects of the finite thickness of turbulent layers of air is of interest where (i) the 
thickness is smaller than the entire atmosphere, cases such as subdivisions in models for multi-conjugated adaptive 
optics or jet-streams, or (ii) where the lateral beam separation is comparatively large, and the scaling down to small 
separations is be done with a single, consistent structure constant, for example comparing interferometric observations 
with atmospheric site monitor data [l^ 

The core of the following is to substitute the approximate calculation of Section 2.2 by transforming (PT|) into (|45l) . 
The nature of the approximation reviewed above is sampling the difference of the two structure functions (jl9p beyond 
the precise limits of the integral; as the structure function is an increasing function of its argument, the net effect of 
the approximation is an exaggeration of "D^p , at least in the case D = Q where the argument of the square root remains 
larger than A/i. 

We first constrain the turbulence model to the limit of infinite outer scale, because the mathematics reduces to 
simpler special functions then (Section IIII Bp . The generic case is phrased in terms of integrals over MacDonald 
Functions in the main part (Section IIII Cp . 



B. In the Kolmogorov Limit 



Taking the limit /o — > 0, (l7|) is inserted into ([T9l) . 

r.i<'/(2sina) /•2h 



_ dAh \^[P^ + {D - AhfY'^ - \Ahpj 

dz £ dy{[l + {D/P~yfy^'^\yr}. 



J -2h 

K Psina rK/iP-H-)) 



sin a K jq 

Binomial expansion of the integrand in hypergeometric series of y transforms 2?^ into 

V^ = ClPP'+-'^g{^,PK) 



(28) 



(29) 



where the dimensionless 



Pi 



K 



P sin a 
K 



is the beam separation in units of the path length, and where 

Pk 



p(i) _ 2^(0) _^ ^(-1) 



with 



(l + 7/2)(l+7)Pi+'' 




-7/2,1/2,1 
3/2,2 




(30) 



(31) 



(32) 
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generalizes (j27l) in terms of three generalized hypergeometric functions. Aspects of the implementation are discussed 
in Appendix \K\ 

Figure [2] illustrates the magnitude of the effect. Plotting 5 as a function of Pk focuses on the relative deviation 
of the full computation in comparison with the standard formula (1261) : non- horizontal lines indicate a nonlinearity 
between I?^ and if/ sin a in (|29p . In general, the phase fluctuations are smaller than the factor 2.9 predicted in the 
approximate ([27|) . which had been derived in the limit K/P — > cx) of infinite layer thickness. At optical and infrared 
wavelengths, the thickness K multiplied by the air mass is limited to the order of 10 km representing the entire 
atmosphere (troposphere), and P is of the order of 2 m for a single telescope up to 100 m for an interferometer. 
Whence Pk is in the range 2 x 10"'' to 0.01, and Figure [2] demonstrates that this finite layer thickness places g in the 
range 2.3 to 2.7. Long baseline interferometry "drops" power of through this geometric sampling effect. 

An increase of D^p is induced if D/P = cotO turns nonzero, i.e., if the pointing direction is not perpendicular to 
the baseline. This tendency is qualitatively understood from Figure [1] since the contribution at D = integrates 
correlations between the two beams at shortest distances P, and nonzero D essentially requires parallel sampling of 
the refractive index structure function correlated over an effective distance b > P. 
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FIG. 2: The factor (|3ip — results in the limit of infinite outer scale — as a function of Pk (the ratio of the lateral distance over 
the path length) for three examples of D/P. 



C. Von-Karman Spectrum 



The intent is to evaluate (ET 



1/Pk 



[P + ifoPn^+'^^' Jo 



dz j dy 



M^TTfy) - Jo(27r/v/l + (D/P - y)2) 



1/Pk 



[/2 + (/oP)2]M/2 7o 

K r pdf 



dz / dy 







2jo(27r/y) - joi2iTf ^iTiD/PWW) 



1/Pk 



dy[l - Pkv] 2joi27r fy) ~ joi2nf^l + [D/P^){i3) 



sin a J, [/2 + (/oP)2](7+3)/2 7„ 
with the convention that ± or ^ indicate summation over both signs in the corresponding term. 



1. First Ray 



The next step is to evaluate the previous equation. It contains a simpler term with a spherical Bessel function with 
argument y, which will be called the contribution of the first ray, and another term depending on D/P, contribution 
of the second ray. Each splits into two terms, one with and one without a factor yPx- Evaluating the /-integrals 
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first, one term is [10|, 3.771.5] [20( 



T3o{27Tfy) 



[/2 + (/oP)2](^+3)/2^ 

f°° tdt 



Its y-integral is 



at 



2^.W^(/opff (3/2 + ,/2) (2^/oi-^)-^-^./.(2^^/oP). (34) 



(Poy)^/'if7/2(^o2/)rfy = ir^-(-Po/-P^)' (35) 

-MJ 



= 7/2, (36) 



where we define [IS 6.561.4] 

K(x)= / t''i^,(t)di = 2^-l^/^^(l/2 + ^y)a;[/^,(a;)L,_l(a;)+/^,_l(a;)L,(x)], (37) 
Jo 

with L representing the modified Struve functions ^, 12.2]. The other contribution of the first ray is P, 11.3.27] 

/ y{P^yr'^K,/^{P^y)dy = --^ [(Po/PK)'+'^i^i+.(Po/P7<) - r(l + v)2^] . (38) 
Jo -fo 

^. Second Ray 

The companion of ([M)) from the second ray is 

[ [/2 + (;ff](..3)/. ^o(2./Vl + bTl)/P]-) 

-[Po Vl + (y T i?/P)2]''i^,(Po + [y T ^/P]2). (39) 



Its y-integral is 



2i+7/2(/pP)7r(3/2 + 7/2)' 



/ (1 - PA-y) [Po ^Jl + iyTD/Pf] ^K, (Po yTTFT^DAPP)dy 

Jo 

= / (1 - PKy T Pa p)[Po Vi + yTK.iPoV^ + y^)dy- (40) 



A substitution Po \/l + = ^ calls to split the interval into the regions y > and y < 0, which is qualified with the 
step-function, <d{x) = 1 if a; > and 8(x) = if x < 0. The contribution from positive y to the integral is 



e(i/PA- T D/p) 



Po 

where 



(1 T pAuAt.P^) + ^r+'K,+,{t) 



p„^i+{i/PKTD/pr 

1 

t^Poy/l+max^{0,TD/P) 



(41) 



C/,(^,w)= / -^L=K,{t)dt (42) 



is an integral described in Appendix |B] The contribution from negative y is 
e(±P>/P) 



Po 



(1 T PK§)C/.(i,Po) - ^t'^+^i^.+i(i) 
P) 



PoVl + (-D/P) = 



t=-PoVl+™2(0,l/PKT-D/P) 



(43) 
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3. Combined Master Equation 



Collecting dSS]), ([38]), (gT]) and (03]) puts ^ into the format 



if 



sma 



5(7,Po,Pif), 



(44) 



with 



(;^)'+'^i^i+.(^o/PK)-r(l + t.)2"^ 



-e{i/PK T ^/P) 



(1 T PK^)U.{t, Po) + ^e+^K,+i{t) 



PoVi + (1/-Pa-T-D/P)2 



-e(±i:'/p) 



(1 T PK§)C/.(t,Po) - ^t-'+'^.+iW 



t=Po^l+max2(0,=F-D/-P) 
PoVl + (^/-P)' 

t=PoVl+min2(0,l/PR-T-D/-P) 



(45) 



at v = 7/2. A summation over the upper and lower sign is still implied for the terms with the two 0-function factors. 
[Note that there are massive cancellations between the group of terms containing and U^, and also between the 
other terms depending on Ki^^, which can be traced back to the differences of the spherical Bessel functions Jq.] 

This is the main result. Equation (|24l) represents the concerted limit Pk — > 0, D/P — > 0, equation ([27]) the limit 
Po ^0, Pk ^ 0, D/P 0, equation ([311) the limit Pq 0, and equation the limit Pq 0, D/P 0. 

Figure [3] shows again g^^'' to point out at which lateral distances the full theory starts to deviate from the straight 
lines predicted by the standard theory. The upper index indicates that this is the contribution from the zeroth order 
of an expansion in powers of (D/P), as in App. |^ 
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FIG. 3: The factor (|45p as a function of the dimensionless Pk for five different dimensionless measures Po of the outer scale, 
with D/P ^0. 



A numerical example: The curve Pq — 0.1 is equivalent to a separation P = 40 cm if the outer scale is I//0 = 25 
m. The standard formula (|24l) assumes a value of 17 = 1.83 at the left edge of the plot; the refined calculation reduces 
this by 25% to g = 1.37 if Psina/K = 0.137, i.e, if the path length through the layer is K/ s'ma 2.9 m. As layers 
are thicker in practise, the more advanced calculation has little impact here. This seems to contradict the importance 
claimed in Section IlII Bl but the following mechanism conciliates both views: As Pq increases, i.e., as the outer scale 
shrinks, the curves in Figure [3] remain horizontal over an increasingly wide initial range on the Pk axis. The standard 
formula remains accurate down to thin turbulent layers, because the clamping of 2?„ introduced by the outer scale 
reduces the importance of correlations over large (including large vertical) scales. The additional strict cut within the 
refined calculation targets smaller values of so the relative effect becomes less important. 
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Inclusion of nonzero "shear" ratios D/P leads to Figure |4l which shows no other features than those expected by 
combining the shapes of Figures [5] and O a hump appears in g if the delay D occupies a major fraction of the path 
length, and g falls off with an approximate — 5/6th power at large Pg- 
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Po=0.001,D/P=0 
Po=0.001,D/P=l 



Po=0.5, D/P=0 
Po=0.5, D/P=l 



FIG. 4: The factor (|45|l as a function of the dimensionless for two different dimensionless measures Po of the outer scale 
and three different ratios D/P. 



IV. SUMMARY 



We developed a formula of the phase structure function depending on the spectral index, on the lateral separation, 
on the air mass and layer thickness, and on the outer scale, which rigorously samples the von-Karman statistics of 
isotropic homogeneous turbulence along two parallel rays. 
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Appendix A: Kolmogorov Limit 

The generalized Hypergeometric Function ([32|) reduces effectively to Gaussian Hypergeometric Functions : 

1/2, 1, -7/2 , \ ^ p / -7/2, 1/2 , \ ^ f -7/2, 1 



( ^' 1/2, 2" " I ^ ) = 2 2^^i ( " U ) - 2^1 ( ^ I z 



3/2 



A useful Taylor expansion of (pij) in powers of _D/P is 
with constant order 

1 



1 + 7/2 



quadratic order 



and biquadratic order 



where 



ff(^) = PK(^-/^-l)(§)^ 



(4) p 7^^/^(2 + 7^-7-V')-^^ iJ^4 



12V'2 
= 1 + 1/Pl^. 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 



Appendix B: Integral U 

If the upper limit z of the integral (|42p is large, it is numerically advantageous to use the complementary [lol . 
6.592.12] 



00 



-.K^{t)dt. 



Further binomial expansion in a power series of w/z. 



leads to a set of integrals where the parameter I is recursively reduced by partial integration if Z > [11 



-K,{z) + 



t-^(i-^)t"-^K,_i{t)dt 



(Bl) 



(B2) 



(B3) 



21-1 ' 1-2Z 

until the case / = is covered by (|37| via [l3, 6.561.16] 

fK^itjdt = y/^2''-^r{u + 1/2) - K(2)- (B4) 

On the other hand, if z is not larger than approximately 1.8w, the repeated partial integration represents a converging 
series 



uaz,u)^j2y 



1 



fc>0 



3-5---(2fc-l)(2fc + l) 



(^2_^2).+ l/2^.-fe^^_^(^)_ 



(B5) 



